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This  paper  develops  an  inverse  problem  approach  to  optimize  the  geometric  structure  of  TECs  (ther¬ 
moelectric  coolers).  The  approach  integrates  a  complete  multi-physics  TEC  model  and  a  simplified 
conjugate-gradient  method.  The  present  TEC  model  couples  the  heat  and  electrical  conductions  and 
accounts  for  all  physical  mechanisms  occurred  within  TECs.  Three  geometric  parameters,  the  semi¬ 
conductor  pair  number,  JV,  leg  length  of  semiconductor  column,  Hpn,  and  base  area  ratio  of  semiconductor 
columns  to  TEC,  y,  are  optimized  simultaneously  at  fixed  current  and  fixed  temperature  difference.  The 
cooling  rate  on  the  cold  end  is  the  objective  function  to  be  maximized  to  obtain  the  optimal  TEC  ge¬ 
ometry.  The  effects  of  applied  current  and  temperature  difference  on  the  optimal  geometry  are  discussed. 
The  results  show  that  at  temperature  difference  of  20  K,  the  geometry  optimization  increases  the  TEC 
cooling  rate  by  1.99-10.21  times  compared  with  the  initial  TEC  geometry,  and  the  optimal  N  decreases 
from  100  to  47  with  invariable  y  =  0.95  and  Hpn  =  0.2  mm,  as  the  applied  current  varies  from  1.0  A  to 
3.0  A.  With  the  increase  in  temperature  difference,  the  optimal  N  increases  at  smaller  currents  of 
/  <  1.0  A,  however,  it  is  almost  invariable  at  larger  currents  of  /  >  1.5  A. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Thermoelectric  devices  can  be  classified  into  two  different 
groups:  one  is  the  TEG  (thermoelectric  generator)  which  converts 
heat  into  electricity  by  Seebeck  effect  [1-3]  and  the  other  is  the  TEC 
(thermoelectric  cooler)  which  converts  electricity  into  heat  by 
Peltier  effect  [4,5].  In  the  recent  years,  as  the  thermoelectric  ma¬ 
terials  have  been  progressed  remarkably,  the  TEC  attracts  more  and 
more  attention  because  it  do  not  use  any  moving  parts  and  envi¬ 
ronmentally  harmful  fluids  as  well  as  it  has  high  reliability  and 
simplicity  [6,7]. 

A  TEC  is  composed  of  a  number  of  semiconductor  element  pairs 
which  are  connected  electrically  in  series  and  thermally  in  parallel, 
and  each  pair  includes  a  p-type  semiconductor  column  and  an 
n-type  semiconductor  column.  The  cooling  performance  of  the  TEC 
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strongly  depends  on  thermoelectric  materials.  The  most  important 
performance  index  of  thermoelectric  materials  is  the  figure  of  merit 
(ZT)  which  is  defined  by  ZT  =  a2Talk ,  where  a  is  the  Seebeck  co¬ 
efficient,  T  is  the  absolute  temperature,  a  is  the  electrical  conduc¬ 
tivity,  and  k  is  the  thermal  conductivity  of  the  material.  A  higher  ZT 
value  will  lead  to  a  better  TEC  performance,  thus,  high  a  and  a,  as 
well  as  low  k  are  desired  for  thermoelectric  materials  [6,7]. 

The  investigations  also  found  that  the  geometric  structure  has 
remarkable  effects  on  the  TEC  and  TEG  (thermoelectric  generator) 
performance  [8-23].  Chen  et  al.  [21]  proposed  a  new  cycle  model 
consisting  of  a  multi-couple  thermoelectric  device  and  involving 
several  key  irreversibilities  of  real  TEGs.  They  used  the  model  to 
optimize  the  performance  of  a  multi-couple  TEG  and  discussed  the 
optimal  structure  of  the  TEG.  Chen  [22]  investigated  the  theoret¬ 
ical  maximum  efficiency  of  solar  TEG  by  using  Lagrangian  multi¬ 
plier  method.  For  a  specific  TEC  with  a  fixed  base  area,  the 
geometric  parameters  include  the  semiconductor  pair  number,  the 
leg  length  (or  height  of  semiconductor  column),  and  the  base  area 
ratio  of  semiconductor  columns  to  TEC.  All  parameters  have 
coupled  effect  on  the  TEC  performance.  Although  an  individual 
parameter  study  is  useful  but  it  cannot  answer  how  one  can  obtain 
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the  optimal  TEC  structure.  Thus,  a  multi-parameter  coupled/ 
combined  effect  should  be  considered  to  reach  the  optimal  TEC 
performance. 

Cheng  et  al.  combined  a  TEC  model  and  a  genetic  algorithm  to 
optimize  the  geometry  and  performance  of  the  single-stage  [23] 
and  two-stage  TEC  [12].  In  their  study,  the  leg  length,  the  leg  area 
and  the  number  of  legs  were  taken  as  the  search  variables  and  were 
optimized  simultaneously  to  reach  the  maximum  cooling  rate  un¬ 
der  the  requirement  of  minimum  COP  (coefficient  of  performance) 
and  the  restriction  on  the  maximum  cost  of  the  material.  In  their 
studies,  the  standard  energy  equilibrium  model  (the  zero¬ 
dimensional  model)  was  used  as  the  direct  problem  solver,  where 
only  the  heat  balance  equations  at  the  hot  and  cold  ends  were 
solved.  The  advantage  of  the  zero-dimensional  model  is  that  it  can 
obtain  the  analytical  expressions  for  the  TEC  performance,  but  its 
accuracy  is  limited  due  to  grossly  simplifying  assumptions.  For 
example,  Joule  heat  and  Thomson  heat  are  assumed  to  be  equally 
distributed  to  the  hot  and  cold  ends,  and  only  the  constant  or 
temperature-averaged  material  properties  can  be  used.  It  is  well 
known  that  the  direct  problem  solver  is  very  critical  for  an  opti¬ 
mization  approach.  If  the  direct  problem  solver  is  not  enough  ac¬ 
curate,  the  optimal  geometry  cannot  be  obtained  by  such 
optimization  approach. 

Recently,  we  developed  a  general,  three-dimensional  numerical 
model  of  thermoelectric  devices  [24],  in  which  the  heat  conduction 
equation  and  the  electric  potential  equation  were  solved  coupled 
and  all  thermoelectric  effects  were  taken  into  account,  including 
Seebeck  effect,  Peltier  effect,  Thomson  effect,  Joule  heating, 
Fourier’s  heat  conduction,  and  heat  loss  to  the  ambient.  The  model 
is  different  from  the  previous  models  [8-20],  in  these  models 
only  the  one-  or  three-dimensional  heat  conduction  equation  with 
Joule  heat  and/or  Thomson  heat  as  internal  heat  sources  was 
solved,  but  electric  potential  equation  was  ignored  based  on  the 
assumption  of  uniform  current  density  across  the  semiconductor 
cross-section. 

The  objective  of  this  work  is  to  develop  an  inverse  problem 
approach  to  look  for  the  optimal  TEC  geometry,  which  integrates 
our  previous  TEC  model  [24]  and  a  simplified  conjugated-gradient 
method.  The  simplified  conjugated-gradient  method  was  originally 
proposed  by  Cheng  et  al.  [25]  and  has  been  applied  successfully  to 
optimize  the  geometric  structures  of  proton  exchange  membrane 
fuel  cells  [26-28],  micro-channel  heat  sinks  [29,30],  and  shell-and- 
tube  heat  exchangers  [31]  in  our  previous  studies.  The  semi¬ 
conductor  pair  number,  the  leg  length,  and  the  base  area  ratio  of 
semiconductor  columns  to  TEC  as  search  variables  are  optimized  to 
reach  the  maximum  TEC  cooling  rate  at  a  fixed  applied  current  and 
a  fixed  temperature  difference  between  the  cold  and  hot  ends.  Then 
the  effects  of  the  applied  current  and  the  temperature  difference  on 
the  optimal  TEC  geometry  are  discussed.  The  present  optimization 
approach  and  the  corresponding  results  are  expected  to  provide  a 
guide  to  the  practical  TEC  design. 

2.  Optimization  approach 

2.1.  Optimized  geometric  parameters 

Fig.  la  shows  the  schematic  of  a  TEC,  which  has  a  base  area  of 
Atotai  and  is  composed  of  N  semiconductor  pairs.  A  single  pair  in¬ 
cludes  a  p-type  semiconductor  column  and  an  n-type  semi¬ 
conductor  column,  three  metallic  interconnectors,  and  two 
electrically  insulating  ceramic  plates.  The  p-type  and  n-type 
semiconductor  columns  have  the  same  leg  length  of  Hpn.  It  is  noted 
that  the  Seebeck  coefficients,  electrical  conductivities  and  thermal 
conductivities  of  p-type  and  n-type  materials  in  the  TEC  are 
different  from  each  other.  Thus,  a  better  TEC  design  should  have 
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Fig.  1.  Schematics  of  TEC  and  single  semiconductor  element  pair:  (a)  TEC;  (b)  single 
semiconductor  element  pair. 


different  cross-sectional  areas  of  p-type  and  n-type  semiconductor 
columns.  However,  the  present  work  does  not  account  for  this 
difference,  and  hence  the  same  cross-sectional  area  of  Ap  =  An  =  Apn 
is  assumed  as  shown  in  Fig.  lb.  The  base  area  ratio  of  semi¬ 
conductor  columns  to  TEC  is  defined  as  y  =  N(AP  +  An)/ 
Atotai  =  (Ap  +  An)/Apair,  where  Apair  is  the  pair  area  and  equals  W0Lo 
(Wo  and  Lo  are  the  width  and  length  of  a  single  pair,  respectively). 
The  copper  interconnector  has  a  thickness  of  Hcu.  Once  the  semi¬ 
conductor  pair  number  N,  the  base  area  ratio  y,  and  the  leg  length 
Hpn  are  given,  the  TEC  geometric  structure  is  determined.  Thus, 
these  three  parameters  are  chosen  as  the  search  variables  and  are 
optimized  simultaneously  to  reach  the  maximum  TEC  cooling  rate. 

2.2.  Direct  problem  model 

The  three-dimensional  TEC  model  proposed  by  our  pervious 
study  [24]  is  adopted  as  direct  problem  model.  The  model  assumes 
that  the  TEC  is  operated  at  steady-state;  the  p-type  and  n-type 
semiconductor  columns  are  two  separate  parts  and  their  material 
properties  are  temperature-dependent;  the  contact  thermal 
resistance  between  semiconductor  columns  and  metallic 
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interconnectors  is  ignored;  and  the  heat  loss  by  radiative  and 
convective  heat  transfer  to  the  ambient  is  ignored.  The  model  in¬ 
cludes  the  heat  conduction  equations  and  the  electric  potential 
equations,  which  are  described  briefly  here  and  more  details  can  be 
found  in  Ref.  [24]. 


K 0  =  ^oTo)  [l  +Ai(T  -  Tq)  +A2(T  -  T0)2]  (7) 

^  =  P(T)  =  P0(T0)[\+B,(T-T0)  +  B2(T-T0)2]  (8) 


2.2.1.  Heat  conduction  equations 


v-(A,vr)+£-ftT-vr  =  o 
v-(Apvr)  +£--/Sp7-vr  =  o 


a(T)  =  ao(T0)[l+Di(T-T0)+D2(T-T0)2]  (9) 

where,  the  parameters  are  fitted  based  on  the  experimental  data  of 

(1 )  Yamashita  et  al.  [32]  and  are  listed  in  Table  1.  The  copper  is  selected 
as  the  interconnector  with  constant  properties  due  to  its  weak 
thermoelectric  effect. 

(2) 

2.3.  Inverse  problem  model 


V-(XnVT) +(*-- PnJ-VT  =  0  (3) 

°n 

2.2.2.  Electric  potential  equation 


V-(<j(V0-aVT))  =  0  (4) 

In  the  above  equations,  A  is  the  thermal  conductivity,  a  is  the 
electric  conductivity,  (3  is  the  Thomson  coefficient,  j  is  the  current 
density  vector,  a  is  the  Seebeck  coefficient,  and  0  is  the  electric 
potential.  The  subscripts  i,  p,  and  n  denote  the  interconnector,  p- 
type  and  n-type  semiconductors,  respectively.  The  Thomson  coef¬ 
ficient  can  be  related  to  the  Seebeck  coefficient  as  follows: 


Once  the  electric  potential  is  obtained,  the  current  density 
vector  can  be  calculated  by  the  following  equation,  or: 

7  =  at  =  o  (-V0  +  aVT)  (6) 

2.2.3.  Boundary  conditions 

The  constant  temperatures  Th  and  Tl  are  applied  at  the  hot  and 
cold  ends  of  the  TEC.  The  temperature  and  heat  flux  are  assumed  to 
be  continuous  on  the  interface  between  the  interconnector  and 
semiconductor.  The  adiabatic  boundary  condition  is  assumed  on 
the  side  surfaces  of  the  TEC.  The  fixed  applied  current  and  zero 
electric  potential  are  specified  on  the  inlet  and  outlet  of  the  TEC, 
and  on  the  other  surfaces  the  current  cannot  flow  out  the  TEC. 

2.2.4.  Material  properties 

The  thermoelectric  materials  Bi2(Teo.94Seo.o6)3  and  (Bio.25Sbo.75) 
Te3  are  chosen  as  the  n-type  and  p-type  semiconductors,  respec¬ 
tively.  Their  thermal  conductivities,  electric  conductivities  and  See¬ 
beck  coefficients  are  all  assumed  to  be  temperature-dependent,  or: 


There  are  two  parameters  used  to  evaluate  the  TEC  perfor¬ 
mance,  the  cooling  rate  Ql  and  the  COP  (coefficient  of  perfor¬ 
mance).  The  cooling  rate  is  defined  as  a  heat  adsorbed  from  the  cold 
end,  and  the  COP  is  defined  as  follows: 

COP  =  ^  ^  (10) 

where  P  is  the  electric  power,  I  is  the  applied  current,  and  V  is  the 
electric  potential  difference  through  the  TEC. 

For  a  TEC  with  a  specific  geometry,  Ql  and  COP  are  all  depen¬ 
dent  on  its  operating  conditions:  the  temperature  difference 
(AT  =  Th-Tl )  and  applied  current.  With  a  fixed  AT,  Ql  and  COP  are 
first  increased  and  then  decreased  as  I  is  increased.  Unfortunately, 
the  optimal  applied  currents  corresponding  to  Ql, max  and  COPmax 
are  not  the  same,  which  means  that  Ql, max  and  COPmax  always 
cannot  reach  simultaneously  [24].  Similarly,  with  the  same  oper¬ 
ating  conditions,  as  the  TEC  geometry  is  varied,  Ql  and  COP  are  all 
varied,  but  maybe  cannot  reach  the  maximums  simultaneously. 
For  practical  cooling  applications,  the  cooling  rate  is  a  priority. 
Hence,  Ql  is  chosen  as  the  objective  function  in  the  present  opti¬ 
mization,  thus,  an  optimal  set  of  geometric  parameters  is  searched 
to  reach  the  maximum  Ql.  Because  optimization  of  TEC  geometry 
maybe  causes  the  reduction  in  the  COP,  the  COP  is  used  as  a 
constraint  condition  during  the  optimization  in  order  to  guarantee 
that  the  TEC  with  the  optimal  geometry  has  a  relatively  high  COP. 
A  value  of  COP  that  is  not  lower  than  70%  COP  of  the  initial  TEC 
geometry  is  chosen  as  the  constraint  condition  here,  or  for  all 
search  steps: 

COP  >  0.7COPini  (11) 

It  is  worth  noting  that  the  selection  of  the  value  of  70%  COPini  seems 
to  be  somewhat  arbitrary,  however,  there  is  no  universal  standard 
to  select  a  proper  COP  when  Ql  is  used  as  the  objective  function 
because  Ql  and  COP  are  two  incomparable  parameters  to  evaluate 
the  TEC  performance. 

The  simplified  conjugate-gradient  method  is  used  as  the  in¬ 
verse  problem  model  in  the  present  optimization  approach.  The 
simplified  conjugate-gradient  method  evaluates  the  gradients  of 


Table  1 

Parameters  in  Eqs.  (7)— (9)  for  thermoelectric  materials  (To  =  300  I<). 


Material 

Thermal  conductivity 

Electric  resistivity 

Seebeck  coefficient 

A0  (W  nr1  K-1) 

A,  (K-1) 

A2  (K-2) 

Po  (□  m) 

B,  (K-') 

b2  (K-2) 

cto  (Vic1) 

D,  (IC1) 

d2  (I<-2) 

p-type 

n-type 

Interconnector 

1.472 

1.643 

400 

-1.29  x  10“3 
-9.80  x  10“4 

0 

1.35  x  10"5 
1.56  x  10"5 

0 

8.826  x  10"6 
8.239  x  10-6 
1.7  x  10"9 

5.88  x  10"3 
4.70  x  10"3 

0 

8.93  x  10"6 
2.67  x  10"6 

0 

2.207  x  10"4 
-2.23  x  lO"4 
6.5  x  10"6 

1.55  x  10"3 
5.62  x  10“4 

0 

-3.15  x 
-4.65  x 

0 

0  0 
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the  objective  function  and  sets  up  a  new  conjugate  direction  for 
the  updated  search  variables  with  the  help  of  a  sensitivity 
analysis.  The  details  of  the  model  can  be  found  in  our  previous 
works  [25-31].  The  optimization  procedure  is  described  briefly 
in  the  following.  For  convenience,  C\,  C2  and  C3  are  used  to 
denote  the  semiconductor  pair  number  N,  the  base  area  ratio  y, 
and  the  leg  length  Hpn,  which  are  referred  to  as  the  search 
variables. 

(1)  Specify  the  initial  guess  for  the  search  variables  C\— C3  and  the 
values  of  the  search  step  sizes  jtfi-fo 

(2)  Create  the  geometry  and  grids  of  the  TEC  using  the  specified 
search  variables  C\— C3.  Set  up  all  boundary  conditions,  and 
then  numerically  solve  Eqs.  (1)— (6). 

(3)  Evaluate  the  objective  function  J  =  Ql-  When  the  convergence 
criterion  is  satisfied,  the  objective  function  reaches  a 
maximum,  then  terminate  iteration;  otherwise,  proceed  to 
step  (4). 

(4)  Evaluate  the  sensitivity  coefficients  6J/0Q  (/<  =  1,  2,  3)  of  the 
objective  function  for  each  search  variable  based  on  the 
following  equation: 


2L 

sck 


AQ 


(12) 


It  is  noted  that  in  step  (4),  the  direct  problem  model  is  invoked 
three  times  to  calculate  (k=  1,2,3),  respectively,  by  introducing 
a  small  perturbation  AQ  into  the  search  variable  Ck  (k  =  1, 2, 3). 

(5)  Evaluate  the  conjugate-gradient  coefficients  yk  {k  =  1,  2,  3)  for 
each  search  variable, 


if - 


-,2 


( dj/dck j*"-1 


(13) 


where,  the  superscripts  n  and  n  -  1  denote  the  nth  and  (n  -  l)th 
search  step,  and  for  n  =  1,  jk  =  0. 


(6)  Evaluate  the  search  directions  izk  (/<  =  1,  2,  3), 


<’  -  (S 


(n) 


(n)  (n-1) 


(14) 


Start 

I 

Initial  guess  parameters 

Q,  Ci  C3 


Grid  generation  unit 

I 

Direct  problem  solver 

I 

Evaluate  objective  function  »  /n) 


Simplified  Conjugate  Gradient  Method 


First  Iteration  yfl)  =0  /=1~3 


Ci=leg  length 

C2=Area  ratio 

C3=Number  of 
semiconductor  pair 

F(C1+AC,,C2,C3)=J i<“> 
Direct  problem  solver 

F(C1,C2+AC2,C3)=^2(n) 
Direct  problem  solver 

F(CuC2,C3+AC$)=J3<a> 
Direct  problem  solver 

|jS 

1 

1 

1 

1 

l 

1 

al 

dCx  AQ 

dC2  A C2 

dC3  A  C3 

Conjugate-gradient 
coefficients  71^ 

Conj  ugate-gr  adient 
coefficients  72(n) 

Conjugate-gradient 
coefficients  7/°) 

n(n)  = 


'a//* /a//°~‘r 

2 

(«)  = 

'a/j'"’  /a/2 

2 

7?- 

dC,  dC, 

•  2 

8C1  dC1 

8C,  8C 

Searching  directions  7Ti(n> 

Searching  directions  x2(n) 

Searching  directions 

(n) 

_(n)  _  __1_  <■)  (n-1) 

1  /1  1 

ocl 

aj  (n) 

^.(n)  _  __2_  +r(n)^(n-D 
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AT  (n) 
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Update  design  variables 
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C,<"*,)=C1<*,-^1<) 

Fig.  2.  Flow  chart  of  optimization  procedure. 
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(7)  Update  new  search  variables, 

c(n+n  =  c(n)  _  pk.n(n)  (15) 

and  then  return  to  step  (2). 

The  typical  values  of  AQ  ( k  =  1,  2, 3)  are  1  for  N,  0.002  for  y,  and 
0.05  mm  for  Hpn.  The  values  of  fa  (/<  =  1,  2,  3)  varies  from  0.001  to 
0.01  depending  on  the  convergence  level  during  optimization.  A 
flow  chart  of  the  optimization  is  shown  in  Fig.  2.  In  steps  (2)  and  (4), 
the  TEC  model,  Eqs.  (l)-(6),  are  converted  to  the  finite-difference 
form  by  the  control  volume  method  and  are  solved  iteratively 
with  an  iteration  criterion  for  convergence  of  10-6.  The  TEC  model 
and  simplified  conjugate-gradient  model  are  coupled  and  solved 
using  a  self-built  Fortran  program. 

The  TEC  model  has  been  validated  in  our  previous  works  for  the 
steady-state  operation  [24]  and  the  transient  operation  [33],  more 
details  can  be  found  in  Refs.  [24,33]. 

3.  Results  and  discussion 

3.1.  The  performance  of  TEC  with  initial  geometry 

A  miniature  TEC  with  the  base  area  of  only  3x3  mm2  is 
selected  as  the  testing  example  to  display  applicability  of  the 
present  inverse  problem  approach  to  optimization  of  multi¬ 
parameters  of  TEC  geometry.  The  initial  design  of  the  TEC  has 
N  =  25,  y  =  0.8,  and  Hpn  =  0.5  mm.  The  cold  end  temperature  TL  is 
fixed  to  300  K.  The  I-Qi ,  /-COP,  and  I—P  curves  for  various  tem¬ 
perature  differences  are  shown  in  Fig.  3.  Fig.  3  confirms  the 
statement  in  Section  2.2,  with  a  fixed  temperature  difference 
AT  >  0,  Qi  and  COP  are  first  increased  and  then  decreased  as  I  is 
increased,  and  there  are  two  optimal  applied  currents  at  which  the 
maximum  Ql  and  the  maximum  COP  occur,  respectively.  To  obtain 
a  higher  COP,  a  small  current  is  recommended,  however,  Ql  is 
lower  (Fig.  3a).  Oppositely,  a  larger  Ql  occurs  at  a  relatively  high 
current,  however,  at  which  COP  is  lower  (Fig.  3b).  Fig.  3  also  shows 
that  Ql  and  COP  at  a  large  AT  are  lower  than  those  at  a  small  AT, 
hence,  the  hot  end  of  the  TEC  generally  is  cooled  by  a  heat  sink  to 
reach  a  better  performance.  The  reduced  COP  at  a  large  AT  can  be 
attributed  to  the  increased  electric  power  P  (Fig.  3c)  as  well  as  the 
reduced  cooling  rate  Q±. 


3.2.  Individual  parametric  study  for  search  variables 


Before  the  optimization,  an  individual  parameter  analysis  is 
performed  to  reveal  the  effects  of  Hpn,  y,  and  N  on  the  TEC  per¬ 
formance,  as  shown  in  Fig.  4.  In  the  calculations,  the  temperatures 
of  the  cold  and  hot  ends  are  all  assumed  to  be  300  K  with  0  K 
temperature  difference.  The  applied  current  is  taken  as  1  A.  The 
parameters  for  the  base-line  case  (corresponding  to  the  initial  TEC 
geometry)  are  Hpn  =  0.5  mm,  y  =  0.8,  and  N  =  25.  In  the  individual 
parameter  study,  when  a  parameter  is  varied  to  evaluate  its  effect, 
other  two  parameters  remain  the  same  as  with  the  base-line  case. 
Fig.  4  shows  that  a  small  Hpn  and  a  large  y  will  elevate  Ql  of  the  TEC, 
however,  an  optimal  N  =  49  is  observed  at  which  the  maximum  Ql 
is  3.52  W.  The  above  results  can  be  explained  using  the  following 
relationship  derived  by  zero-dimensional  model  [34], 


Ql=N 


_Hpn\ 

G  n^pn  / 


( ^p^pn  AnApn\ 

V  Hpn  Hpn  ; 

(16) 


where,  the  first,  second,  and  third  terms  of  the  right  side  in  Eq.  (16) 
denote  the  heat  absorbed  from  the  cold  end  due  to  Peltier  effect, 


Fig.  3.  Performance  of  TEC  with  initial  geometry:  (a)  cooling  rate;  (b)  COP;  (c)  electric 
power. 


the  heat  transferred  to  the  cold  end  due  to  Joule  heating,  and  the 
heat  transferred  to  the  cold  end  due  to  Fourier’s  backward  heat 
conduction.  According  to  Eq.  (16),  when  the  temperature  difference 
is  AT  =  0  K,  the  third  term  equals  zero,  hence,  the  cooling  rate  Ql  is 
determined  by  subtracting  the  second  term  from  the  first  term, 
thus,  the  increase  in  the  first  term  and/or  the  decrease  in  the  second 
term  will  elevate  Q±.  With  the  same  N  and  the  same  operating 
conditions  (/  =  1  A,  TL  =  300  K,  and  AT  =  0  K),  a  small  Hpn  and  a 
larger  y  (corresponding  to  a  large  cross-sectional  area  of  semi¬ 
conductor  column  Apn)  mean  a  low  electric  resistance  of  the  TEC 
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Fig.  4.  The  effect  of  individual  parameter  on  TEC  performance:  (a)  leg  length  Hpn;  (b) 
base  area  ratio  y;  (c)  semiconductor  pair  number  N. 


and  hence  a  low  electric  potential  difference  through  the  TEC 
(Fig.  5),  thus,  the  Joule  heating  is  weak  which  leads  to  a  large  Q±. 
However,  with  the  same  Hpn  and  y,  a  large  N  corresponds  to  a  small 
y4pn,  hence,  the  Peltier  heat  N(ap-an)ITi  and  the  heat  transferred  to 
the  cold  end  0.5  NHpn  l((JpApn)+'l l{(TnApn))  are  all  increased  as  N 
is  increased.  Therefore,  increase  in  N  has  both  positive  and  negative 
effect  on  the  cooling  rate,  in  other  word,  the  cooling  rate  is  non¬ 
monotonic  function  of  N,  thus,  the  optimal  N  corresponding  to 
the  maximum  Ql  is  dependent  on  competition  of  the  first  term  and 
the  second  term.  Different  from  Ql,  COP  is  always  varied  monoto¬ 
nously,  a  higher  COP  is  observed  at  smaller  Hpn  and  N,  as  well  as  at 
larger  y. 


3.3.  Three-parameter  optimization  of  TEC  geometry 

The  individual  parameter  study  indicates  that  a  small  Hpn  and  a 
large  y  improve  Q±  of  the  TEC,  however,  Hpn  cannot  be  arbitrarily 
small  due  to  the  constraint  by  practical  processing  technology,  and 
y  cannot  be  equal  to  1,  hence,  a  constraint  condition  with 
Hpn  >  0.2  mm  and  y  <  0.95  is  adopted  in  the  present  optimization. 
It  can  be  expected  that  the  optimal  geometry  of  the  TEC  should  be 
varied  at  different  operating  conditions  due  to  the  dependence  of 
Ql  and  COP  on  the  applied  current  and  temperature  difference. 
Therefore,  the  geometric  parameters  (N,  y,  Hpn)  of  the  TEC  are 


Fig.  5.  Electric  potential  difference  through  the  TEC  for  various  geometric  parameters: 
(a)  leg  length  Hpn;  (b)  base  area  ratio  y;  (c)  semiconductor  pair  number  N. 
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optimized  to  look  for  the  maximum  Ql  at  operating  conditions  of 
I  =  0.5  A,  1.0  A,  1.5  A,  2.0  A,  2.5  A,  and  3.0  A  as  well  as  AT  =  0  K,  20  K, 
40  K,  and  60  K,  respectively.  It  is  worth  noting  that  although 
AT  =  0  K  is  no  practical  meanings  for  refrigeration  purpose,  how¬ 
ever,  it  can  obtain  the  maximum  Ql  for  cooling  purpose  at  the  same 
TEC  geometry  and  the  same  applied  current. 

3.3.1.  Optimal  design  at  AT  =  0  I< 

Fig.  6  shows  the  changes  in  three  search  variables  for  various  /  at 
AT  =  0  K  during  optimization.  It  can  be  seen  that  the  search  vari¬ 
ables  vary  from  their  initial  values  of  Hpn  =  0.5  mm,  y  =  0.8,  and 
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Fig.  6.  Variations  of  search  variables  during  optimization:  (a)  leg  length  Hpn;  (b)  base 
area  ratio  y;  (c)  semiconductor  pair  number  N. 


N  =  25  as  the  search  step  number  increases,  and  finally  reach  the 
optimal  values.  Optimal  design  has  the  same  y  =  0.95  and 
Hpn  =  0.2  mm,  however,  has  different  N  at  various  applied  currents 
/.  It  is  noted  that  y  and  Hpn  for  optimal  design  are  equal  to  their 
constraint  value,  hence,  when  the  base  area  of  the  TEC  is  fixed,  to 
obtain  the  optimal  cooling  rate  the  leg  length  should  be  as  small  as 
possible  and  the  total  area  of  semiconductor  columns  should  be  as 
large  as  possible  for  a  practical  TEC  design.  Fig.  6c  shows  that  the 
optimal  N  is  100  at  /  =  0.5  A,  and  it  reduces  to  95,  90,  72,  58,  and  48 
at  I  =  1.0  A,  1.5  A,  2.0  A,  2.5  A,  and  3.0  A,  respectively.  This  result  can 
be  explained  as  follows.  Since  the  optimal  design  has  the  same 
y  =  0.95  and  Hpn  =  0.2  mm,  a  larger  N  means  that  the  semi¬ 
conductor  column  has  a  smaller  cross-sectional  area  Apn,  which 
causes  a  higher  electric  resistance  of  the  TEC  and  hence  a  stronger 
Joule  heating.  This  effect  becomes  more  significantly  at  higher  /,  the 
Joule  heating  reduces  Ql  based  on  Eq.  (16),  thus,  the  optimal  design 
at  higher  /  approaches  to  reduce  N  to  obtain  the  maximum  Ql,  as 
compared  to  at  lower  /. 

The  variations  of  Ql  and  COP  for  various  I  during  optimization 
are  shown  in  Fig.  7.  Ql  increases  with  the  search  steps,  it  approaches 
a  plateau  of  5.56  W  after  15  search  steps  for  /  =  0.5  A,  about  276% 
higher  than  that  for  the  initial  design.  For  currents  I  =  1.0  A,  1.5  A, 
2.0  A,  2.5  A,  and  3.0  A,  the  optimal  design  has  Ql  =  9.12  W,  10.36  W, 
10.47  W,  10.46  W,  and  10.46  W,  about  252%,  214%,  198%,  243%,  and 
624%  higher  than  that  for  the  initial  design,  respectively. 

The  present  optimization  restricts  COP  at  each  search  step  to  be 
not  lower  than  70%  COP  of  the  initial  design  (Eq.  (11 )).  Fig.  7b  shows 
that  COP  of  the  optimal  design  at  small  currents  of  /  <  1.5  A 
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Table  2 

The  optimal  semiconductor  pair  numbers  for  different  temperature  differences  and 
applied  currents. 


/(A) 

0.5 

1.0 

1.5 

2.0 

2.5 

3.0 

Ar  =  0  K 

100 

95 

90 

72 

58 

48 

AT  =  20  K 

- 

100 

88 

71 

57 

47 

AT  =  40  K 

- 

110 

89 

70 

56 

46 

AT  =  60  K 

- 

110 

90 

69 

56 

- 

approaches  this  limit  value.  However,  at  large  currents  of  I  >  2.0  A, 
the  optimal  design  has  not  only  higher  Qi,  but  also  higher  COP  than 
the  initial  design.  It  is  noted  that  there  is  a  peak  value  of  COP  during 
optimization  as  shown  in  Fig.  7b,  which  means  that  even  for  small 
applied  currents,  COP  and  Ql  can  also  be  simultaneously  higher 
than  those  of  the  initial  design.  For  example,  at  I  =  0.5  A,  the  initial 
design  has  COP  =  3.96  and  Ql  =  1.48  W,  the  design  at  tenth  search 
step  has  COP  =  5.25  and  Ql  =  3.28  W,  and  the  design  at  twelfth 
search  step  has  COP  =  4.13  and  Ql  =  4.21  W.  These  two  designs  have 
lower  Ql  but  higher  COP  than  the  optimal  design,  hence,  if  large 
COP  and  Ql  are  needed  to  be  considered  simultaneously  for  a 
specific  application,  these  design  can  be  selected. 

3.3.2.  Effect  of  temperature  difference  on  the  optimal  design 

Apart  from  the  applied  current  /,  the  optimal  TEC  design  is  also 
dependent  on  the  temperature  difference  AT.  The  three-parameter 
optimizations  are  also  carried  out  at  AT  =  20  K,  40  K,  and  60  K  with 
a  fixed  cold  end  temperature  of  TL  =  300  K.  Similarly,  the  optimal 
design  has  Hpn  =  0.2  mm  and  y  =  0.95.  However,  the  optimal 
semiconductor  pair  number  N  is  different  for  various  AT  and  /,  as 
shown  in  Table  2.  At  the  same  AT,  a  large  N  occurs  at  low  applied 
current,  and  it  decreases  with  the  current  increased,  this  phe¬ 
nomena  has  been  explained  in  Section  3.3.1.  At  I  <  1.0  A,  a  large  N 
occurs  at  high  temperature  difference,  however,  at  I  >  1.5  A,  N  is 
almost  independent  on  the  temperature  difference.  This  result  is 
particularly  useful  for  practical  applications,  since  a  sole  optimal 
design  can  obtain  the  maximum  cooling  rate  at  various  tempera¬ 
ture  differences. 

Figs.  8  and  9  compare  COP  and  Ql  of  the  optimal  and  initial 
designs  at  AT  =  0  K,  20  K,  40  K,  and  60  K.  Again,  COP  of  the  optimal 
design  is  lower  at  small  currents  I  <  1.5  A,  but  is  higher  at  large 
currents  I  >  2.0  A,  as  compared  to  that  of  the  initial  design  (Fig.  8); 
However,  Q±  is  significantly  increased  for  all  temperature  differ¬ 
ences  and  applied  currents  (Fig.  9).  It  is  noted  that  at  the  same  AT, 
when  I  >  2.0  A  almost  the  same  COP  and  Q±  are  obtained  for  the 


optimal  design.  Though  the  present  optimization  is  carried  out  for  a 
TEC  with  a  specific  base  area  of  3  x  3  mm2,  but  the  results  obtained 
here  provide  a  design  direction,  which  is  useful  to  guide  the  prac¬ 
tical  TEC  design. 

4.  Conclusions 

This  work  develops  an  optimization  approach,  which  combines 
a  simplified  conjugate-gradient  method  and  a  complete  three- 
dimensional  TEC  model,  to  look  for  optimal  geometric  structure 
of  TECs.  A  miniature  TEC  with  a  base  area  of  3  x  3  mm2  is  used  to 
test  the  optimization  approach.  The  semiconductor  pair  number  N, 
the  leg  length  Hpn,  and  the  base  area  ratio  y  are  adopted  as  search 
variables  to  maximum  the  objective  function,  defined  as  the  cool¬ 
ing  rate  of  the  TEC.  A  value  of  COP  that  is  not  lower  than  70%  COP  of 
the  initial  design  at  each  search  step  is  used  as  a  constraint  con¬ 
dition  during  optimization.  The  initial  design  has  N  =  25, 
Hpn  =  0.5  mm,  and  y  =  0.8.  The  optimizations  are  carried  out  in  a 
wide  range  of  operating  conditions:  at  temperature  difference 
AT  =  0  K,  20  K,  40  K,  and  60  K,  as  well  as  at  applied  current  I  =  0.5  A, 
1.0  A,  1.5  A,  2.0  A,  2.5  A,  and  3.0  A. 

The  optimal  design  has  the  same  Hpn  =  0.2  mm,  and  y  =  0.95  for 
all  AT  and  /,  which  reach  their  limit  values  constrained  by  the 
present  work,  hence,  when  the  base  area  of  the  TEC  is  fixed,  to 
obtain  the  maximum  cooling  rate  the  leg  length  should  be  as  small 
as  possible  and  the  total  area  of  semiconductor  columns  should  be 
as  large  as  possible.  The  present  optimizations  also  show  that  a 
large  N  is  better  to  increase  the  cooling  rate  of  the  TEC  at  small 
currents  and  high  temperature  differences.  At  small  currents  of 
I  <  1.5  A,  improved  cooling  rate  for  the  optimal  design  is  at  the 
expense  of  reduced  COP,  however,  at  large  currents  of  /  >  2.0  A,  both 
the  cooling  rate  and  the  COP  can  be  increased  significantly  by 
searching  for  the  optimal  design. 

Although  the  combined  optimization  is  proven  effective  only  for 
the  miniature  TEC  with  the  base  area  of  3  x  3  mm2,  it  is  expected 
that  the  proposed  design  strategy  is  likely  applicable  to  designs  of 
the  practical  TEC  with  larger  base  area. 
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